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I.  INTRODUCTION 


The  continuing  need  for  experimental  damage  and  effects  data  for 
high  energy  laser  interactions  has  resulted  in  a number  of  experimental 
programs  within  the  three  services.  The  outputs  of  these  programs  are 
compilations  of  burn-through  or  failure  times  as  a function  of  such 
parameters  as  target  material,  target  surface,  shield  thickness,  air 
flow,  and  particularly  beam  size,  shape  and  intensity.  Unfortunately, 
the  available  laser  devices  which  are  powerful  enough  for  realistic 
experiments  permit  only  very  imperfect  control  of  beam  size  and  average 
intensity,  with  no  control  over  shape  and  the  intensity  of  spikes.  As 
a result,  experimentalists  are  forced  to  set  up  their  tests  based  on 
an  anticipated  average  intensity,  carefully  diagnose  the  beam  actually 
delivered,  and  interpret  their  results  after  the  fact. 

For  continuous  wave  lasers,  the  best  available  technique  for  beam 
diagnostics  involves  the  use  of  an  infra-red  (IR)  scanning  camera  in 
conjunction  with  a real  time  beam  power  measurement  (e.g.  a chopper) 
and  ball-calorimetry  for  absolute  calibration.  The  scanning  camera 
measures  the  diffuse  reflection  of  the  beam  from  a specular  surface  by 
imaging  the  (diffuse)  beam  reflection  and  sequentially  scanning  the 
image  in  a line-by-line  fashion,  similar  to  the  scan  of  a black-and- 
white  television  camera.  If  a turning  flat  is  used  as  the  specular 
surface,  the  specularly  reflected  beam  can  be  directed  into  the  test 
target  while  the  diffuse  image  is  being  IR-scanned;  thus,  a real-time 
measurement  is  made.  The  IR-scan  results  are  only  relative,  since  the 
shape  of  the  beam  is  shown,  but  not  the  absolute  intensities.  However, 
absolute  values  are  easily  derived  by  integrating  over  the  entire 
scanner  picture,  arriving  at  the  scanner-relative  total  power,  and 
normalizing  to  the  real-time  absolute  total  power  measurement. 

The  results  obtained  from  this  normalization  are  absolute  two- 
dimensional  maps  of  intensities  through  the  beam  profile  at  the 
specular  surface.  Unfortunately,  the  profiles  encountered  are  usually 
so  erratic  that  the  mere  reporting  of  one  test,  let  alone  test-to-test 
correlation,  is  extremely  cumbersome.  Therefore  a need  exists  to 
identify  an  "ideal-beam' equivalent"  for  the  actual  beam  of  each  test. 

The  physical  output  of  an  IR-scanning  camera  is  an  analog 
(magnetic)  tape.  This  memorandum  will  trace  the  techniques  and  soft- 
ware developed  by  BRL  to  take  the  analog  information  and  arrive  at 
contour  plots,  "equivalent  beams",  and  three-dimensional  prospective 
pictures . 
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II.  DATA  REDUCTION  PROCEDURE 


The  steps  in  the  data  reduction  sequence  are  shown  in  the  flow 
diagram  in  Figure  1. 

In  step  1,  the  analog  tape  is  digitized.  The  analog  tape  is  format- 
ted as  follows:  Each  frame  begins  with  a frame  synch  pulse;  similarly 

each  line  scan  (22  lines  on  the  BRL  IR  Scanning  Camera)  begins  with  a 
line  synch  pulse.  These  synch  pulses  are  superimposed  on  the  continuous 
analog  signal,  with  a frame  synch  pulse  occuring  every  0.008  seconds  for 
the  BRL  scanner.  Each  line  is  then  sampled  at  42  equally  spaced 
intervals  in  an  apparent  22  x 42  matrix.  However,  timing  problems  result 
in  the  first  line  being  entirely  dropped  (much  like  the  retrace  blanking 
is  done  on  a television  receiver).  Further,  a line  synch  pulse  ringing 
problem  necessitates  zeroing  out  the  first  three  samples  on  each  line. 

The  resulting  scanner  matrix  is  thus  21  x 42,  with  non-zero  data  in 
21  x 39  samples. 

The  digitizing  for  BRL-processed-scanner  data  is  accomplished  in 
a mobile  data  acquisition  system.  The  digitizing  procedure  results  in 
a 9-track  digital  tape  consisting  of  seven  packed  16-bit  words  followed 
by  882  8-bit  bytes.  The  16-bit  words  contain  such  information  as  the 
run  and  frame  number.  The  digitizing  procedure  also  produces  identifi- 
cation of  the  frame  numbers  corresponding  to  the  start  and  stop  of  the 
laser  beam. 

Step  2 consists  of  taking  the  9-track  digitized  data  tape  and 
transferring  the  information  to  a 7-track  tape  which  can  be  accommodated 
on  a UNIVAC  1108  computer.  This  is  merely  a system  requirement  and  is 
accomplished  by  established  system  routines;  no  data  reduction  is 
involved. 

In  step  3 the  program  VLC  performs  the  normalization  discussed  in 
Section  I,  and  produces  a line-printer  contour  plot  of  the  normalized 
data.  The  time  integrated  raw  data,  plus  certain  terms  are  also  dumped 
onto  a binary  tape.  In  step  4,  the  binary  tape  conveys  the  necessary 
data  to  the  smoothing  and  fitting  program,  DRAPER. FIT.  The  results  of 
DRAPER. FIT  are  available  in  synopsized  form  from  the  line  printer.  A 
comparison  of  the  smoothed  data  to  the  DRAPER. FIT  Gaussian  is  available 
in  3-dimensional  prospective  CALCOMP  pictures  from  the  program  3D-PLOT, 
shown  as  step  5.  The  VLC,  DRAPER. FIT,  and  3D- PLOT  programs  are  discussed 
in  the  succeeding  sections. 


III.  PROGRAM  VLC 

The  program  VLC  performs  five  functions:  (1)  reads  the  digitized 

data  tape  and  converts  it  to  computer  manageable  integers;  (2)  determines 
the  background  noise  and  subtracts  it  from  the  raw  data;  (3)  normalizes 
the  raw  data,  as  discussed  in  section  I;  (4)  outputs  line-printer  contour 
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plots;  and  (5)  produces  a binary  tape  for  use  in  step  4,  Figure  1.  The 
program  consists  of  two  routines:  the  main  routine  and  IRTAPE,  an 
assembly  language  tape  interpreter.1 

The  inputs  to  VLC  are  shown  in  Table  1.  The  data  tape  must  be 
called  INPUT,  the  print  file  is  unit  10,  and  the  binary  output  tape  is 
unit  12. 

Referring  to  Table  1,  the  TEST  Ident  is  any  verbal  phrase  to  appear 
as  an  output  heading.  FOV  (field  of  view  of  the  scanner)  is  the  linear 
width  (=  height)  of  the  scanner  picture  on  the  turning  flat.  Since  the 
FOV  is  normally  unchanged  for  several  runs,  the  last  inclusive  run 
number  for  which  a given  FOV  is  to  apply  is  input  along  with  the  FOV. 

A new  FOV  requires  a new  card  input  but  not  until  after  the  data  for  the 
previous  FOV's  runs  has  been  input. 

BEAM  ON  and  OFF  FRAMES  are  obtained  with  the  digitized  tape  and 
input  as  shown  on  Figure  1.  The  beam  on/off  card  is  also  used  to  submit 
a start  frame  if,  for  example,  the  early  frames  are  destroyed.  For 
such  cases,  a $ in  column  1 signifies  that  3 integers  follow. 

The  "number  of  time  increments  and  values"  card  serves  one  of  two 
functions;  it  either  indicates  the  number  of  sequential  time  integrals 
to  be  made  and  their  time  intervals,  or  signals  a calorimeter  run.  The 
latter  requires  1,  0.  on  the  time  card;  the  time  integral  is  then  taken 
over  the  entire  set  of  beam  on/off  frames.  If  the  time  card  is  1,  X 
where  X does  not  equal  0.,  a component  run  is  assumed.  In  either  case, 
if  only  one  time  is  given,  no  further  cards  are  required  for  that  run. 

On  the  other  hand,  if  more  than  one  time  is  given,  the  next  card 
identifies  the  shot  sequence;  e.g.  shield,  component  #1,  component  #2. 
(This  information  is  only  used  for  headings  and  does  not  affect  the 
data  reduction.) 

A detailed  key  to  the  hard  copy  output  of  the  VLC  program  is  given 
in  Appendix  A. 


IV.  SMOOTHING  AND  FITTING  THE  MEASURED  BEAM 


A.  SMOOTHING 

As  discussed  in  section  I,  the  beams  being  diagnosed  are  often 
pathologically  indescribable.  The  irregular,  unsymmetrical  shapes  are 
dotted  by  spikes  containing  little  power  but  having  peak  intensities 
many  times  higher  than  the  local  average.  An  example  of  such  a beam 
is  shown  in  figure  2A  and  2B.  From  an  effects/damage  experimental 
viewpoint,  however,  these  spikes  - if  real  - are  often  irrelevant. 


lMr.  John  Kinah , USABRL } private  communication . 
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FIGURE  1.  FLOW  DIAGRAM 


10 


Table  1.  VLC  Inputs 

Card  No. 

1.  @XQT  AMUCK  * SCANNER. VLC 

2.  BRL  SERIES  3A  - TEST  Ident. 

3.  4.75,  2180  - FOV,  last  inclusive  RUN  NO. 

4.  2175  - RUN  NUMBER 

5.  33,  296*  - BEAM  ON  and  OFF  FRAMES 

6.  75.3  - TOTAL  POWER 

7.  2,  0.5,  1.2  - Number  of  time  incr  § values 

8.  2,  3 - Code  corresponding  to  times 

(USED  ONLY  IF  .GT.  1 time  is  given  above) 
where  1 = CALORIMETER 

2 = SHIELD 

3 = COMPONENT  #1 

4 = COMPONENT  #2 

For  a calorimeter  test  in  which  the  beam  intensity  is  to  be  averaged 
over  the  entire  shot,  the  time  card  (#7)  should  be  1,0.  since  the  number 
of  time  increments  is  only  1 and  the  code  card  (#8)  is  omitted. 

9.  2176  - Same  as  #4 

10.  29,  421  - Same  as  #5 


29.  2180  - Same  as  #4 

30.  37,  226  - Same  as  #5 

31.  81.2  - Same  as  #6 

32.  1,  0.5  - Same  as  #7 

33.  5.25,  2191  - Same  as  #3 

@E0F 

UNIT  CONTAINS 

INPUT  Data  tape 

10  Print  file 

12  Output  binary 

* Optional  input  for  START/STOP  card  is  $,  33,  296,  10  - BEAM  ON  and  OFF, 
START  FRAME 
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Figure  2A.  Measured  Beam  Observed  from  the  '21'  Direction 

Perspective  of  square  f ield-of-view  displayed  21  lines  x 42  samples. 
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(Note:  Unpredictable  errors,  such  as  small  irregularities  on  the  mirror 

surface,  could  cause  apparent  spikes.  Discussion  of  experimental  error 
is,  however,  beyond  the  scope  of  this  report.)  Certain  phenomena  tend 
to  spread  out  the  beam  power  upon  interaction  with  the  target.  Since 
the  spikes  contain  little  power,  their  effect,  once  spread,  is  minimal. 

It  is  therefore  appropriate,  as  part  of  scanner  data  reduction,  to  show 
the  effect  of  this  spreading  by  applying  a suitable  smoothing  to  the 
raw  data. 

Often  the  most  significant  spreading  phenomena  is  radial  conduction. 
Therefore,  to  provide  a lower  bound  on  the  smoothing,  a simple  radial 
conduction  model  was  constructed.  The  model,  presented  in  Appendix  B, 
indicated  that,  for  a typical  experiment,  resolution  below  0.5  cm.  at 
the  target  is  probably  superfluous.  However,  the  model  and  attendant 
computer  code  can  be  run  on  a case  by  case  basis. 

The  DRAPER. FIT  program  uses  this  resolution  information.  The 
analyst  is  queried  for  "X-cells  over  which  to  smooth"  (recall  - X is 
the  line  L direction;  i.e.,  the  21  direction  in  the  21  x 42  output 
matrix.)  The  number  of  cells  is  found  by  determining  a minimum 
resolution  from  the  conduction  model,  and  dividing  that  by  the  cell 
size  as  output  by  VLC.  DRAPER. FIT  then  replaces  each  cell  value  by  a 
weighted  2-dimensional  average  of  local  cells.  The  effect  of  averaging 
the  beam  of  Figure  2 over  one  and  two  cells  is  shown  in  Figures  3 and  4 
respectively. 

B.  FITTING 

As  seen  from  the  above  figures,  the  resulting  beam  is  still  not 
conveniently  described.  Test-to-test  comparison  of  results  is 
impossible  because  of  the  inability  to  quantitatively  account  for  the 
important  beam  parameters.  What  is  needed  is  an  equivalent  beam  of  a 
standard  form  to  be  equated  to  each  measured  actual  beam.  For  operational 
reasons,  a Gaussian  is  a preferrable  standard  form.  The  task  then 
reduces  to  defining  a technique  for  fitting  an  equivalent  Gaussian  to  the 
smoothed  beam  profile  data;  i.e.,  finding  the  appropriate  parameters  for 
the  form: 


I(x,y)  = Io*exp(-(x-xo)2/2Rx2)*exp(-(Y-Yo)2/2Ry2)  + B (1) 


The  crux  of  the  problem  lies  in  determining  what  is  meant  by 
"equivalent".  For  instance,  a CHI-square  minimization2  program  was  used 
to  fit  equation  1 to  some  actual  data.  It  was  found  that  the  area  away 

2J.T.  Klopaia , "FNFIT:  An  Easy-To-Use,  Arbitvary  Funation-To-Data  Fitting 
Routine  (A  User's  Manual),"  U.S.  Army  Ballistic  Research  Laboratory , 

BRL  MR  2402  (1974).  (AD  #A000655) 
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from  the  peak  dominated  the  fit,  resulting  in  values  for  I which  were 

clearly  too  small  to  be  indicative  of  the  actual  beam.  The  I was 

improved  by  weighting  each  data  point  by  Z(I)a,  where  Z(I)  is  the 
(intensity)  value  at  the  point,  and  a is  a positive  number.  As 
expected,  the  value  of  I more  closely  approached  the  height  of  the 

dominant  peaks  as  larger  values  for  a (e.g.,  1.0,  1.5,  2.0,  3.0)  were 
investigated.  Unfortunately,  the  arbitrary  choice  of  values  for  a was 
unacceptable,  and  other  definitions  of  equivalent  were  sought. 

Another  possible  approach  to  defining  an  equivalent  Gaussian  is  to 
determine  various  characteristics  of  the  actual  beam,  and  choose  a 
Gaussian  having  the  same  values  for  those  characteristics.  Examples  of 
such  characteristics  are  the  mean  and  standard  deviation.  If  as  many 
(independent)  characteristics  are  chosen  as  there  are  parameters  in  the 
Gaussian  (six  in  equation  1),  then  the  Gaussian  is  uniquely  determined. 

A possible  set  of  six  characteristics  is  shown  in  Table  2. 


Table  2.  Beam  Characteristics  (Power) 


1.  j I(x,y,t)-dA-dt  = C • j dA-dt,  A'  outside  2ox,  2 ay 
A'  A' 


2,  3.  ux 

4,  5.  a2 
a2 

y 

6.  E 


/ 

/ 

/ 

/ 


X' I (x,y,t)dA*dt/E;  = / y- I (x,y,t) -dA-dt 

(x-yx)2  • I(x,y,t) -dA-dt  / E 
(y-Py)2  * I(x,y,t)*dA*dt  / E 
I (x,y,t) *dA*dt 


Characteristic  1 is  a definition  of  background,  and  establishes  a 
baseline  level  (C) . Characteristics  2 and  3 define  the  mean,  and  4 and 
5 define  the  variance.  E is  recognized  as  the  total  energy  delivered. 
Therefore,  characteristic  6 assures  that  the  total  power  of  the  equiva- 
lent Gaussian  is  equal  to  the  (time  averaged)  total  power  of  the  beam. 
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The  integrals  involved  above  can  be  done  in  closed  form  for  the 
Gaussian  of  equation  1,  yielding  the  expressions  for  the  six  character- 
istics shown  in  Table  3. 


Table  3.  Gaussian  Values  for  Beam  Characteristics  (Power) 
1.  B = C 

2,  x0 » Y0  = means  (yx,u  ) 

4,  5.  R2,  R2  = variances  ( o2, 

x y x 

6.  2*ir*R  »R  • I • t = E 
x y o 


The  POWER  option  of  the  DRAPER. FIT  program  calculates  the  equivalent 

Gaussian  via  these  six  criteria.  The  characteristics  C,  y , y , a , a , 

x y x y 

and  E are  evaluated  from  the  digitized  frames  by  performing  the  needed 
spatial  and  temporal  integrations  numerically.  The  relations  shown  in 
Table  3 are  then  used  to  compute  the  equivalent  Gaussian  parameters 


R 


I and  B. 
o 


The  results  of  fitting  the  beams  of  Figures  2,  3 and  4 via  this 
technique  are  shown  in  Figures  5,  6,  and  7.  DRAPER. FIT  also  gives 
information  printed  along  the  left  margin  of  the  picture  (see  Appendix  C) . 
The  order  and  meaning  of  the  information  is  given  in  Appendix  C. 

Finally,  DRAPER. FIT  also  gives  a CHI-SQ  for  each  fit,  where  CHI-SQ  is 
defined  as 


NDATA 

CHI-SQ  = N * I (Z  - I(X  , Y ))2/I2(X  Y.)  (2) 

1=1 

where  N = 1 / (NDATA  - 7) 


For  the  fits  in  Figures  5,  6 and  7,  the  CHI-SQ  values  were  0.644, 
0.471  and  0.404  respectively.  No  absolute  interpretation  for  CHI-SQ  is 
put  forth  here  in  terms  of  statistical  confidence  levels;  rather, 

CHI-SQ  is  used  as  a relative  measurement.  (Recall  that  the  CHI-SQ 
minimization  approach  was  unsatisfactory.)  As  such,  the  CHI-SQ  indicates 
the  improvement  in  fit  with  smoothing,  as  seen  in  Figures  5,  6 and  7. 
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Figure  5A.  DRAPER. FIT  (Power  Option)  for  Figure  2A 

Perspective  of  square  field-of-view  displayed  21  lines  x 42  samples. 
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Figure  6A.  DRAPER. FIT  (Power  Option)  for  Figure  3A 

Perspective  of  square  field-of-view  displayed  21  lines  x 42  samples 
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Figure  6B.  DRAPER. FIT  (Power  Option)  for  Figure  3B 

Perspective  of  square  field-of-view  displayed  21  lines  x 42  samples. 
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Figure  7A.  DRAPER. FIT  (Power  Option)  for  Figure  4A 

Perspective  of  square  field-of-view  displayed  21  lines  x 42  samples. 
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The  choice  of  characteristics  upon  which  to  base  the  fit  is,  as 
implied  above,  somewhat  arbitrary.  For  the  reporting  of  damage  and 
effects  data,  the  "DRAPER  characteristics",  a set  slightly  different 
from  that  in  Table  2,  is  preferred. 

In  a report  on  parameter  sensitivity  analysis  considerations  by 
the  Charles  Stark  Draper  Laboratories,3  a discussion  of  beam  character- 
ization was  presented.  In  particular,  it  was  noted  that  the  term  I(x,y,t) 
serves  as  a distribution  function  for  the  calculation  of  x and  s (y  and  a). 
One  is  then  led  to  define  an  average  intensity,  I,  analogously. 

I = j I (x,y,t) • I (x,y,t) *dA*dt  / E (3) 

The  Draper  report  went  on  to  discuss  the  behaviour  of  this  I for  a set 
of  idealized  beams,  showing  that  the  damaging  potential  of  the  beams  was 
appropriately  ranked  by  I. 

The  Draper  observations  were  embraced  by  the  equivalent  Gaussian 
portion  of  the  scanner  program.  I constitutes  a characteristic  of  the 
beam  in  the  spnse  of  those  in  Table  2,  with  the  added  benefit  that  I 
is  a reflection  of  damage  potential.  Furthermore,  substituting  I for  E 
(characteristic  6 in  Table  2)  results  in  a complete  and  mutually 
independent  set,  as  listed  in  Table  4. 


Table  4.  Beam  Characteristics  (Draper) 


1. 


j I (x,y,t) -dA*dt  = C • j dA*dt,  A' 


outside  2o  , 2o 
x’  y 


2,3.  u 


_ / x-I(x,y,t)dA-dt/E;  vy  - / y* I (x,y,t) *dA*dt 

4,5.  a2  = J (x-Mx)2  • I (x,y,t) *dA*dt  / E 

°y  = f ('y’yy')2  ’ 1 ’dA’dt  / E 

6.  I = J [I (x,y,t)] 2*dA*dt  / E 


3 Olssonj  Ouellette t et  al.3  "Parameter  Sensitivity  Analysis 3"  the 
Charles  Stark  Draper  Laboratory 3 Ina.3  SEE # 009-05-S46  (June  1975), 
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As  in  the  POWER  normalized  case,  the  expression  for  I for  the 
Gaussian  of  equation  1 is 

! = f-rr-R  *R  *I2*t)  / E = I / 2 (4) 

v x y o } o' 

The  counterpart  to  Table  3 is  presented  as  Table  5. 


Table  5.  Gaussian  Values  for  Beam  Characteristics  (DRAPER) 


1.  B ~ C 


2,3.  X , Y = means  (u  , u ) 
o o x y 

4,5.  R2,  R2  = variances  (a2, a2) 

x y xs  y 

6.  1=2*1 
o 


The  beams  of  Figures  2-4  are  presented  with  their  DRAPER. FIT  in 
Figures  8-10.  As  expected,  the  DRAPER. FIT  results  more  closely  fit 
the  peaks.  The  overall  fits,  as  determined  by  the  CHI-SQ  criterion, 
are  better,  especially  for  the  unsmoothed  beam.  CHI-SQ  values  for 
Figures  8,  9 and  10  are  0.558,  0.449,  and  0.394  respectively. 

The  DRAPER. FIT  computer  program  is  an  interactive  routine  which 
solicits  desired  inputs  in  free  field  format.  The  input  routine  directly 
reads  the  output  tape  from  VLC  program.  Outputs  of  the  DRAPER. FIT 
program  include  line  printer  copy  (alternate  print  file  7 optional),  and 
a tape  compatible  with  the  3D-PL0T  routine. 


V.  PROGRAM  3D- PLOT 

The  3D-PL0T  routine  does  no  data  processing.  However,  as  can  be 
seen  by  comparing  Figures  2-10  with  Figure  A-l  in  Appendix  A,  the 
interpretation  of  scanner  data  is  markedly  assisted  by  graphical  displays. 
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ure  9A.  DRAPER. FIT  (Gaussian)  for  Figure  3A 

Perspective  of  square  field-of-view  displayed  21  lines  x 42  samples 


c/) 

<D 


CM 

Tf 


X 


{/l 

<D 

C 


CQ 

O* 

0) 

fH 

3 

bfl 

•rH 

Uh 


36 


37 


Figure  9C.  DRAPER. FIT  (Gaussian)  for  Figure  3A  Oriented  Clockwise  45  Degrees 

Perspective  of  square  field-of-view  displayed  21  lines  x 42  samples. 
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iRAPER.FIT  (Gaussian)  for  Figure  4A 

erspective  of  square  field-of-view  displayed  21  lines  x 42  samples. 
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pective  of  square  field-of-view  displayed  21  lines 
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3D-PL0T  is  based  on  a routine  originally  prepared  by  the  National 
Center  for  Atmospheric  Research4;  however,  extensive  modifications  were 
made  to  allow  proper  handling  of  the  sharp  peaks  and  irregular  projection 
angles  which  can  occur  in  scanner  data.  The  data  points  (X,  Y,  Z,  F, 
where  F is  the  fit)  are  read  in  any  order  from  a tape  via  NTRAN  commands, 
which  allows  tape  positioning  under  program  control.  The  left  margin 
information  is  conveyed  from  the  DRAPER. FIT  program  via  the  data  tape. 

If  both  the  data  and  the  fit  are  to  be  plotted,  two  colors  are  used 
(This  two  color  contrast,  unfortunately,  does  not  show  in  Figures  2-10.) 

All  inputs  are  solicited.  These  include  input  and  output  units, 
scale  factor,  eye  position,  focal  point  and  tape  positioning  instructions. 
The  scale  factor  is  used  to  convert  the  Z and  F values  (still  in  counts  - 
see  explanation  of  PEAK  = in  Appendix  A)  to  a visually  pleasing  size  with 
respect  to  the  X and  Y dimensions.  To  aid  in  the  determination  of  the 
scale  factor,  the  PEAK  = counts  to  PK100  kW  factor  (PKCON)  is  output  by 
the  DRAPER. FIT  program.  Hence,  using  PKCON  as  the  scale  factor  will 
result  in  the  peaks  being  one  unit  high  in  the  Z direction  for  every 
kW/cm2,  as  compared  to  X and  Y of  21  and  42,  respectively.  Scalefactor  = 
10*PKCON  was  used  in  Figures  2-10. 


VI . SUMMARY 

The  hardware,  software  and  experience  exists  to  process,  interpret 
and  display  IR-scanner  data.  The  codes  are  production  oriented, 
resulting  in  fast  turnaround  time  with  a minimum  of  tape  handling.  The 
slowest  portions  of  the  chain  shown  in  Figure  1 are  the  physical  transport 
ing  of  the  tape  to  and  from  step  1 (step  1 is  presently  accomplished  at 
a site  remote  from  BRL)  and  the  delays  involved  with  producing  CALCOMP 
drawings . 


4Afr.  Thomas  Wright , National  Center  for  Atmospheric  Research } Boulder } 
Colorado. 


41 


ACKNOWLEDGEMENT 


The  authors  wish  to  acknowledge  the  help  of  Mr.  Earl  Weaver,  the 
originator  of  VLC,  and  especially  Mr.  John  Kinch,  who  arranged  steps 
1,  2 and  3 in  the  program  flow  (Figure  1).  Mr.  Kinch  also  wrote  the 
machine  language  programs  to  unfold  the  digitized  tapes,  and  rewrote 
VLC,  making  it  Univac  compatible  and  extending  its  capabilities. 


42 


APPENDIX  A 

POWER  DISTRIBUTION  CONTOUR  PLOT 


The  computer  printout  in  Figure  A-l  shows  a composite  frame  of  beam 
diagnostic  data  for  a particular  test  plus  relevant  frame  parameters. 

The  composite  frame  is  generated  by  summing  over  the  test  time  the 
power  density  in  each  cell  of  each  frame  comprising  the  composite, 
normalizing  to  the  (input)  total  power,  and  then  printed  as  an  array  of 
numbers  ranging  between  0 and  9 to  indicate  the  contour  level  of 
accumulated  power  in  each  cell.  The  following  descriptions  apply  to  the 
hard  copy  output  of  the  VLC  program. 

PROJECT:  computer  storage  and  user  identity  code 

TEST  NO.:  device  run  number 

TYPE:  target  classification;  i.e.,  calorimeter  or  shield  or  component 

1 or  component  2 

TIME  X.XXX:  test  time  in  seconds  for  the  target  TYPE  (determined  from 

bum  wires  or  high  speed  photographic  measurements).  Times 
are  input  into  the  computer  code. 

XXX  FRAMES  FROM  XXX  ON:  contour  plot  data  is  a composite  representing 

XXX  frames  of  scanner  data  beginning  with  frame  XXX.  The 
'start'  frame  number  (also  the  'stop'  frame  number)  is 
obtained  from  the  analog  tape  record.  The  start  and  stop 
frame  numbers  are  inputs  to  the  computer  code. 

NOISE  IS  X COUNTS:  Number  of  noise  counts  in  each  cell  of  those  frames 

used  to  generate  the  composite  frame.  Noise  count  is 
determined  for  each  cell  from  the  frames  preceeding  the 
'start'  frame  and  substracted  from  each  cell  to  remove  noise 
level  from  the  contour  data.  If  noise  counts  are  divided 
by  the  number  of  counts  in  PEAK  = XXXX,  signal-to-noise 
ratio  for  peak  power  can  be  estimated. 

Count  is  defined  as  the  (arbitrary)  number  of  bits  of 
information  assigned  to  the  voltage  value  of  the  peak  power 
by  the  digitizing  circuitry  and  is  the  reference  for  scaling 
all  power  amplitude  data. 

TOTAL  POWER:  total  test  power  in  kW  (as  measured  by  the  chopper-detector 

beam  sampler  of  the  device) . TOTAL  POWER  is  an  input  to 
the  computer  code. 

L AXIS:  line  axis  of  21  rows 

S AXIS:  sample  axis  of  42  columns 
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***  u.  s.  army  ballistic  RESEARCH  LABORATORIES  ** 
VULNERABILITY  and  LETHALITY  DIV.  rad.  eng.  br 
data  REDUCTION  DATE:  122776  1610^1 


va- 
e chac 

O CC  IT  I 
• • 

\L 

c 


ID 

©c  r-  r*  o 

c r c i 


K. 

c 


cr*- 
r-  o <r  ^ o 
hie  I 

• i |T 

— c\ 


E 

nl 

<u 

CQ 


K)  O O 
o cr  CV  O 

u c + 
• • r- 

— < lO 


^4  1/5 


* o 

IT  O ff  H C 
lOtfK)  + 


Cl 

E 

O 

u 


z 

o 


c 
« c 


y 

c cr 
^ cr  uj 
iru  1 
h-  o 
Z in  Q. 
LlJ  UJ 

Z * -J 
o < < 
c.  cr  k 
y u o 
o ►— 

O IT 
OJ 


l/) 


UJ 

Q. 


1/5 

K 

z 

Z5 

o 

u 


-J 

a 


OJ 

< 

© 

X 

UJ 

OJ 

• 

• • 

c 

h- 

IT 

CJ 

UJ 

»- 

“> 

1/1 

c 

UJ 

cr 

>— 

a 

4-4  LJ 

h*  tn 


ccc 
ccc 
O C o 
o c © 
o c c 
oco 
o o o 
oc  c 
o c e 
c c o 
o o o 
o c c 
© O •-* 
O C «H 
O O *“♦ 
© © ^ 
O o *H 

o 

H o ^ W 

O OJ 
x O O *h 

o c c 
«c  o w ** 

O o -H 

© c o 
o o o 
o o o 
o o o 
oco 
o o o 
o © c 
o o o 
o o o 
o o c 
© o © 
o o o 
o o o 
o cr  o 
o o o 
o o o 
o o o 
o o o 


ccc 
o o o 
oco 
o e o 
oco 
o o c 
oco 
oco 
coo 
occ 
o o o 
O O «H 
HO  (\i 
C «H  Oi 

O O «H 

—4  o o 


HH  OHM 


HK)  H 

W ^ C\J 

H(MlO 

hH  OJ 


in 


^ OJ  lO 
-*  K)  lO 
*H  — 4 OJ 
—4  OJ 
H lO  H 
HIOH 
O — 4 — < 
o H H 
O *H  •“! 
O O «H 
O O «H 
O © «H 

o o o 
o o o 
o o o 
o o o 
o o o 
o o o 


ccc 
oco 
o o o 
o c c 
ccc 
o o o 
O O o 
o o c 
o o c 
o c c 

© C — 4 
•-» 

M M H 

OJ  r 
CVJ  Ol  CM 
HMitf 
H (\l  ^ 

ro  IT  OJ 
lO  IT  lO 
OJ  OJ  ^ 
OJ  \C  ^ 
(VI  (T  ^ 

H ^ <C 
W + 
— 4 OJ  ID 
OJ  ^ ID 

iom^ 

— 4 OJ  O) 


c c 
o o 
o o 
c o 
c c 
o o 
o o 
o o 
c c 
o c 
o o 

C «H 
O 


H (VI  H 


•H  O 
O *H  o 

O^O 
© — 4 O 

o o o 
o o o 
oco 
o o o 
o o o 


—4  fO 
OJ  OJ 
OJ  — < 
OJ  *H 

K)  0.1 

if  ^ 

oj  r» 
K)  in 

O'  OJ 
« lO 
© Ki 
CG  OJ 
fO  IO 
»H  to 
—4  OJ 


4 OJ 


—4  O 
—4  O 

o o 
o o 
o o 
o o 
o o 
o o 


c c c c 
o o o o 
o c o o 
o c o o 
c c o c 
o o o o 
o o o o 
o o o o 
o c o o 
c o o o 
o o o o 
o c o o 
*H  o o o 
^ *H  O o 

•H  O © O 
—4—4  0 0 

W «H  o o 

NHHO 
OJ  w w o 
(V)  HH  O 
OJHH  O 
'H  H H O 
—4000 

W -4  O O 
w4  -4  O O 
W *H  O O 
N H O O 

ki  -<  o © 

MHO  O 
-4000 
—4000 
—4000 
—4000 

c o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 
o o o o 


c c c c 
o o o c 
o o o c 
o o c c 
c c c c 
o o o o 
o o c o 
o o c o 
o o c c 
CO  c c 
o o o o 
©ccc 
o o o o 
c o o c 
o o o o 
© o o o 
o o o o 
© o o o 
o o o o 
o o o o 
o © o o 
coco 
o o o o 
© o o o 
© o © o 
© o o o 
© © o o 
O © © c 
o o o o 
o © © o 
© o o o 
© © © o 
© o o o 
o o o o 
o o o o 
© o o o 
© o o o 
o o o © 
o o o o 
o o o o 
o o o o 
o o o o 


O ID  o 

? C CC  K C 
if  ID  C ♦ 

• • H 

If  o 


Ol  ^ C 

m o O'  oj  c 

K)MH  + 

• • KT 

o o- 

—4  vD 


vC 

•3- 


< 

UJ 

0. 

II 


o 

cr 


u 

II 

♦ 


S'  • 

c 

II 

OJ  < 
H UJ 

• cr 

OJ  < 


in 

a • 

-I 


oj  c a 
O •- 


J c if 

in  in 

»o  ^ 


Jh 

o 

4-1 

3 

a 

E 

o 

u 


i*  O 

r-  o 
c + 
• © 
o © 
OJ  o 


Off  H 

— 4 in  -3-  o 
© —•  + 
• • o 

a.  c 


«C  X H tsi 


II 

o 

-J 

J 

— . — 

OJ 

vO  Cf 

o 

—1  *H 

O O' 

rO 

o O O ID 

• 

l/l  <C  CO 

® OJ  Ou 

c cr 

io  r- 

o o-  in 

• • 

II 

CO  O' 

OJ  O'  O' 

K 1 

i/i 

• • 

• • • 

ID 

in 

OJ 

•“* 

® 

UJ 

in 

ISJ 

1 HH 

H 

l/l 

a:  x:  uj 

D 

O o 11 

D < -J  S’ 

I 

a x k 

a o ® x 

K UJ  UJ  ^ 

J _J 

o < z 

o *h  o <c 

za  > uj  j 

i UJ 

O UJ  UJ  > 1C  1C  u 

o uj  cr  o 

u 

U CL  O 

< CL  0,  G. 

un  J<> 

44 


+ = CENTROID:  location  of  the  cell  containing  the  center  of  mass 

of  the  total  power 

* = PEAK:  location  of  the  cell  containing  the  highest  accumulation 

of  power  density  (always  a '9'  cell) 

FIELD:  size  of  the  field  of  view  of  the  scanner  shown  in  the  composite 

frame . 

CELL  SIZE:  'S'  and  ’ L ' dimensions  of  resolution  cell  in  centimeters  and 

area  in  cm2 

COORD:  coordinates  of  PEAK  and  CENT  given  in  linear  dimensions  and 

column/row  designation  with  respect  to  the  origin  positioned 
in  the  upper  left-hand  corner  of  the  array 

AVGPD:  average  power  density  of  all  cells  in  the  array  with  numbers 

greater  than  zero;  i.e.,  the  TOTAL  POWER  divided  by  AREA  of  the 
cells  of  CONTUR  1 

PK100:  100  percent  of  the  peak  power  density  defined  as  the  accumulated 

power  density  in  the  PEAK  cell  (kW  per  cm2) 

PK98 : 98  percent  of  the  peak  power  density 

PEAK  = XXXX:  number  of  counts  in  the  cell  containing  peak  power  density 

CONTUR  0,1, 2,.... 9:  contour  intervals  corresponding  to  the  contour 

numbers  in  the  array 

% PEAK:  minimum  value  of  the  percentage  of  peak  power  density  contained 

in  the  contour  interval;  e.g.,  for  contour  2,  all  cells  of  at 
least  20%  of  the  peak  power  density  but  less  than  30%  are 
designated  as  2 in  the  array 

LEVEL:  value  of  the  power  density  at  the  beginning  of  the  contour 

interval  designated  by  % PEAK;  e.g.,  for  CONTUR  4 (40%)  and 

40 

PK  100  = 5.00,  LEVEL  would  be  5.00  x = 2.00 

AREA:  area  of  that  portion  of  the  array  within  the  contour  designated 

by  CONTUR  as  well  as  all  higher  numbered  contours 

VOLUME:  total  power  in  kW  within  AREA 
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APPENDIX  B 

RADIAL  CONDUCTION  LIMITS 


SPECIAL  RESOLUTION  REQUIREMENTS  FOR  HEL  BURN  DIAGNOSTICS 
We  assume  that  damage  is  thermal . 

Rationale:  Consider  two  thermal  spikes  located  a distance  2X 

apart,  each  containing  energy  Q. 

|«-  X -*•  «-  X +| 


If  one  wishes  to  resolve  these  spikes,  the  diagnostic  tool  would 
require  spacing  between  resolution  elements  of  .LT.  X.  The  question 
is  how  small  is  X still  meaningful?  Our  answer:  X is  certainly 

meaningless  for  values  smaller  than  Xq,  where  XQ  is  the  minimum  resolution 
limit  imposed  by  other  causes. 

A major  contributor  to  Xq  is  dictated  by  radial  conduction.  A heat 
conduction  equation  was  employed  to  determine  the  temperature  difference 
between  T and  T (i.e.  (T  - T )/T  , precisely),  as  a function  of  X, 

O X 0X0 

time,  spike  to  ambient  ratio,  and  thinness  and  material  of  plate  being 
irradiated.  For  these  calculations,  the  following  assumptions  were 
made: 


(1)  the  problem  could  be  treated  as  1-D.  This  corresponds 
to  spikes  that  have  a large  Y dimension  relative  to  X. 

(2)  the  plate  is  thermally  thin  to  the  beam,  and  the  power 

is  therefore  distributed  throughout  the  thinness.  (NOTE:  The  error  in 

this  assumption  is  largely  mitigated  by  the  fact  that  ratios  of  differ- 
ences are  being  taken.)  The  power  in  the  spike  per  unit  Y is  C ^ gp j " 

Iamb) * whereas  ambient  power  per  unit  Y is  Hence  the  ratio 

of  the  power  in  the  spike  to  ambient  is 

('spike  - W*THIN''IAMB*X  3 ( RATIO- 1)*THIN/X 

Here,  THIN  is  the  thickness  of  the  plate. 

Outputs  from  the  computer  code  are  shown  in  Figures  B-l  through 
B-8.  The  units  are  centimeters  and  seconds,  and  the  diffusivity  is 
appropriate  for  aluminum.  As  expected,  the  importance  of  the  spike 
decreases  with  increasing  time,  and  increases  with  increasing  X,  RATIO, 
and  plate  thickness.  Notice,  too,  that  RATIO  and  THIN  are  not  independent 
the  curves  are  actually  a function  of  (RATIO  - 1)  * THIN/2.  Hence 
Figure  B-7  is  equivalent  to  Figure  B-4. 
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It  remains  to  apply  these  results  to  the  question  of  beam 
diagnostics.  A 0.127  cm  (50  mil)  sheet  is  an  average  shield.  Consider 
an  average  flux  of  3 kW/cm2,  with  spikes  as  high  as  9 kW/cm2.  Burn 
through  should  take  in  the  order  of  0.10  second.  Error  in  the  time 
measurement  is  probably  ^20%.  Figure  B-8  gives  a 20%  temperature 
difference  at  0.4  cm.  Hence,  it  makes  no  sense  to  sample  more  closely 
than  that.  A more  significant  difference  occurs  at  0.6  cm.  At  0.5  cm 
for  X,  the  ratio  of  power  in  the  spikes  to  ambient  is  0.5,  which  is 
quite  pathelogical  (hefty  spikes!).  It  therefore  appears  that  0.5  cm 
is  a lower  limit  on  the  reasonable  resolution  at  the  target. 

The  following  descriptions  apply  to  the  hard  copy  outputs  of  the 
computer  code  for  temperature  difference  in  Figures  B-l  through  B-8: 

DIFFUSIVITY  = X.XXX  rate  of  diffusion  of  aluminum  plate 

SPIKE/AMBIENT  = X.XX  ratio  of  power  in  the  spike  to 

ambient  power  per  unit  Y 

THINNESS  OF  PLATE  = X.XXX  thickness  of  plate  (aluminum)  in  cm 
TEMP  (0,T)  temperature  at  spike  at  time  T 

TEMP  (L,  T)  temperature  between  spikes  at  Time  T 

TIME  .050  to  .500  sec  at  dt  = .050  sec 

DISTANCE  TO  BOUNDARY  .200  to  2.000  cm  at  0.2  cm  intervals 
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APPENDIX  C 
DRAPER. FIT 


The  composite  frame  data,  as  shown  on  the  contour  plots  (reference 
Appendix  A),  is  smoothed  to  reflect  the  resolution  limits  imposed  by 
radial  conduction.  It  is  then  fit  by  an  equivalent  Gaussian  of  the  form 

10, y)  = IQ*exp  (-(x-xo)2/2o2)  *exp(-(y-y0)2/2cr2)+C 


Here  x and  y are  the  energy  distribution  centroid,  o and  a the 
o Jo  x y 

standard  deviation  of  the  distribution,  and  I is  the  Equivalent 

Gaussian  intensity.  C is  the  baseline.  This  intensity  is  derived 
either  from  an  I2  integration  based  on  an  intensity  criterion  appearing 
in  a report  by  the  Charles  Stark  Draper  Company  or  on  total  power.  The 
following  descriptions  apply  to  the  hard  copy  output  of  the  DRAPER. FIT 
program  shown  in  Figure  C-l. 


***RUN  NUMBER  XXXX:  device  run  number 

TOTAL  POWER  XXX. XX:  total  test  power  in  kW  (measured  by  the  chopper- 

detector  beam  sampler  of  the  device) 


NOISE  X:  number  noise  counts  per  cell 

X TIME  INTERVALS. . .X. XXX  X.XXX:  no.  of  time  intervals  for  the  run 

number  with  duration  of  time  intervals  in  seconds 


FIT  # XX:  sequence  number  assigned  by  the  Draper  program  for  either 

POWER  or  DRAPER .FIT 

CHI  SQ.  = X.XXX:  results  of  chi-square  test  fitting  the  composite 

frame  contour  data  against  the  Gaussian  distribution 
by  DRAPER  program 

PKCON  = X.XXX-XX:  conversion  factor  is  ratio  of  PK100  to  PEAK  = used 

to  convert  the  Z scale  of  the  3D-PL0T  program  to 
compatible  range  with  respect  to  the  S and  L axes. 

(The  following  information  also  appears  on  3D-PL0T  Calcomp  plots) 

RUN  NO.  XXXX.X:  run  number  and  time  sequence  number  assigned  to  FIT  # 

above 


ON  TIME  X.XX:  test  time  in  seconds 

TP  XXX. X:  total  test  power  in  kW 
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Figure  Cl.  Computer  Printout  - DRAPER. FIT  Program 


N X:  number  of  noise  counts  in  RUN  NO.  above 


(XX. X,  XX. X):  L and  S coordinates  of  the  peak  value  of  the  Gaussian  fit 

X.X.X.X:  standard  deviation  of  Gaussian  fit  in  L and  S directions  in  cm 

X.XXX:  peak  value  of  the  Gaussian  fit  in  kW  cm  2 

.XXX:  baseline  of  Gaussian  fit  in  kWcm-2 
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ATTN:  Mr.  G.  McLafferty 

Mr.  A.  Angelbeek 
400  Main  Street 
East  Hartford,  CT  06108 

1 Vought  Corporation 

ATTN:  Mr.  B.  Forcht 

P.  0.  Box  5907 
Dallas,  TX  75222 

2 Battelle  Memorial  Institute 
Columbus  Laboratory 

ATTN:  Mr.  K.  Wilkes 

Mr.  S.  Rubin 
505  King  Avenue 
Columbus,  OH  43201 

4 Massachusetts  Institute  of 
Technology,  Lincoln  Lab 
ATTN:  Dr.  Marque t 

Dr.  Rediker 
Dr.  Edelberg 
Dr.  Reis 
P.  0.  Box  73 
Lexington,  MA  02173 


No.  of 

Copies  Organization 

1 Stanford  Research  Institute 
ATTN:  Dr.  R.  Armistead 

Menlo  Park,  CA  94025 

1 Stanford  Research  Institute 
ATTN:  Mr.  John  H.  Hennings 

306  Wynn  Drive,  NW 
Huntsville,  AL  35805 


Aberdeen  Proving  Ground 
Dir,  USAMSAA 

ATTN:  DRXSY,  Dr.  Sperrazza 

DRXSY-AFF , Mr.  D. Smith 
Mr.  A.  Henderson 
Marine  Corps  Ln  Ofc 
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